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Abstract 

In an extreme mass-ratio binary black hole system, a non-equatorial orbit will list (i.e. increase its 
angle of inclination, i) as it evolves in Kerr spacetime. The abutment, a set of evolving, near-polar, 
retrograde orbits, for which the instantaneous Carter constant (Q) is at its maximum value (Qx) 
for given values of latus rectum (I) and eccentricity (e), has been introduced as a laboratory in 
which the consistency of dQ/dt with corresponding evolution equations for dl/dt and de/dt might 
be tested independently of a specific radiation back-reaction model. To demonstrate the use of the 
abutment as such a laboratory, a derivation of dQ/dt, based only on published formulae for dl/dt 
and de/dt, was performed for elliptical orbits on the abutment. The resulting expression for dQ/dt 
matched the published result to the second order in e. We believe the abutment is a potentially 
useful tool for improving the accuracy of evolution equations to higher orders of e and l" 1 . 
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I. INTRODUCTION 



An extreme mass-ratio binary black hole system (EMRI) is composed of a primary object, 
which can be a Kerr black hole of mass M ~ 10 6 — 10 7 solar masses with a spin 1 S = |J| /M 2 
(where J is the spin angular momentum), and an orbiting secondary object of mass m ~ 
1 — 10 solar masses. Theoretical models to describe the orbital evolution of the secondary 
object in various situations have been derived and presented in the literature: circular orbits 
in the equatorial plane of the primary object [THE], elliptical orbits in the equatorial plane 
[71-113]. and an extensive body of research on circular or elliptical orbits inclined with respect 
to the equatorial plane [T4T131j. Such models are used to generate hypothetical gravitational 
waveforms (GW), which provide templates for use in the detection of gravitation wave 
signals by pattern recognition (Punturo et al. [32J). The detection of GW radiation by 
the Earth-based Laser Interferometer Gravitational Wave Observatory (LIGO) or the Laser 
Interferometer Space Antenna (LISA) depends fundamentally on the availability of correct 
templates [I3J E3J E3 • 

Performing direct observations of relativistic effects is an important challenge. The So- 
lar System affords one the opportunity to observe and model the motions of natural and 
artificial bodies in Kerr spacetime in the weak-field, slow motion limit [35, 36J; and recent 
measurements of artificial-satellite orbits have produced estimates of the Lense-Thirring 
precession to an accuracy of 10% [35J. Further, the discovery of Sagittarius A*, a massive 
black hole (MBH) of ~ 4.0 x 10 6 solar masses, at the centre of our galaxy (see [37-39] and 
references therein), offers a new opportunity to study Kerr spacetime by the observation of 
various stars in inclined, highly elliptical orbits, and by the analysis of their orbital dynamics 
[37M40] . Relativistic effects are difficult to discern since the orbital periods of the stars are 
in the tens of years |40j . and for orbits that come close to the MBH, tidal disruption is a 
concern ([10] and see Appendix B in [37]); yet, observation has great potential to aid in 
the study of Kerr spacetime. In the case of an EMRI, unfortunately, the part played by 
a theoretician is a fiduciary one; thus, the introduction of tools with which the evolution 
equations can be tested for consistency is most beneficial: the abutment is one such tool, 
but it is not intended to replace existing methods. 

1 Our use of S for black hole spin arose during our initial studies of the work of Barack and Cutler |24| . 
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The concept of the abutment, a boundary that defines a set of near-polar retrograde 
orbits, was developed and introduced by P. G. Komorowski in his Doctoral thesis [41 J and 



in a previous work [42j (we shall review the abutment in detail in section II B 1); two uses 
of the abutment had emerged: first, it suggested a means of testing the consistency of the 
evolution of the Carter constant of circular orbits (dQ/dt) with respect to that of the latus 
rectum (dl/dt); and second, it permitted a numerical analysis of the rate of change of the 
orbital angle of inclination, l, with respect to I (idi/dl) ) for circular orbits constrained 

V / min 

to evolve along the abutment. In this work we shall extend these uses to orbits of non-zero 
eccentricity (0 < e < 1) by testing the consistency of expressions for dQ/dt with expressions 
for dl/dt and de/dt, and we shall perform an analytical treatment of i and the list rate of 
the same. Further, a physically realistic orbital evolution follows the abutment (Qx) in only 
one case, the evolution of an orbit in a Schwarzschild black hole (SBH) system (S = 0). We 
shall now consider the general case of an evolving orbit that intersects the abutment, Qx, 
tangentially at a single point (contact of the first order (see 99 in |43j)) as it follows a path 
defined by Q pa th- Further, by performing our analysis for elliptical orbits, the abutment 
becomes a two dimensional surface that defines the maximum value of Q for given values of 
e and latus rectum, I = l/M . Therefore one must view the abutment as a set of contiguous 
points rather than a path to be followed by an evolving orbit; and it is at these points that 
the derivatives, dQx/dl and dQx/de, fix the corresponding slopes of Q pa th- But as reported 
in |42j . the second-order effect 2 must be included when working with l at the abutment. 

In section [TT] we shall analytically derive the formula for i for elliptical orbits on the 
abutment, and thus confirm the result for idc/dl J |42j . which was derived numerically 

V / min 

for circular orbits. In addition, we shall analytically derive di/de for elliptical orbits that 
evolve on the abutment. In section [TTT1 we shall include the effect of the second derivative 
°f Qpath (i- e - the second-order effect) by introducing reductive ansatze for circular and 
elliptical orbits, and thus create a more physically realistic model for an evolving orbit at 
the abutment. 

Because our abutment model is independent of any specific radiation back-reaction model, 
we now have a laboratory that allows us to perform tests of established listing formulae. In 



2 When we refer to the second-order effect at the abutment, we refer to the second derivative of Qx with 
respect to I, not to S 2 . See [H] for background discussion. 
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section |IV[ we shall demonstrate the usefulness of the abutment in testing the consistency 
of dQ/dt equations with respect to dl/dt and de/dt evolution equations, and in calculating 
di/dt for elliptical orbits of small eccentricity (i.e. near-circular). In section |V| we shall 
conclude our work and recommend directions that warrant further study. 

We define i to be the maximum polar angle reached by the secondary object in its 
orbit (see equation (42) in [42J). This definition differs from that used by others (Gair and 
Glampedakis [27] and Glampedakis, Hughes, and Kennefick [23J); but when performing our 
analysis to the leading order in S, there is no significant difference. 

II. AN ANALYTICAL FORMULA FOR THE ANGLE OF INCLINATION OF AN 
ELLIPTICAL ORBIT ON THE ABUTMENT 

A. Introduction 

The listing of an inclined elliptical orbit of eccentricity (e) can be described by dijdl and 
dt/de, where i is the angle of inclination of the orbit and I is its latus rectum normalised 
with respect to the mass (M) of the Kerr black hole (KBH). A set of essential analytical 
formulae for the orbital constants of motion has been derived in [42] : the Carter constant at 
the abutment (Qx), the orbital energy (E), and the quantity, X = L z — SE, as well as an 
analytical formula for l in terms of these constants of motion. Numerical analysis yielded 
an equation for ydt/dl\ for circular orbits: 

/ min 

- (l22.75 - 36S 3 ) r 11/2 - (63/25 + 35/4S 3 ) T 9/2 

min 

- 15/2St 7/2 - 9/2St 5/2 . (1) 

To verify equation ([!]) analytically, we shall derive the result to order 3 in S (i.e. 0(S 3 )). 
Observe that equation ([!]) is a series expansion in terms of f~5. Further, the series coefficients 
are themselves series expansions of odd powers of S. These are important properties, which 
we shall confirm and investigate. Equation Q is not sufficient for understanding the effect 
of radiation back-reaction on the listing of near-polar orbits; therefore, it is necessary to 
develop an analytical formula for 1 on the abutment so that a more thorough treatment 
can be made. We shall review the analytical formulae reported in |42| for elliptical orbits, 
and develop appropriate expansions of those formulae in terms of S. The MacLaurin series 
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expansions of the functions 1/(1 + x), y/1 + x, arccos(x), cos(x), and sin(x) are essential for 
this work. 



B. Review of Analytical Formulae 

1. The abutment, Qx 

The analytical formula for X± (where X = L z — SE) for elliptical and inclined orbits 
about a KBH was found to be [42]: 

V 2 Z$ + Z^Q ± 25 V Z 7 ZgZg 

A ± ~ ~rz ~ 2 \ ' W 

!(l[3-l + e 2 ) -4S 2 (l-e 2 ) 2 ) 

where 

Z 5 = P { (l + 3e 2 + l) S 2 - l(3 - 1+ e 2 ) } , (3) 
Z & = -2 (1 - e 2 ) 2 5 4 + 21 (2e 4 + (2 - l\ e 2 + 4 - f) S 2 - I 2 (3 - 1 + e 2 )' , (4) 
Z 7 = S 2 (l + e) 2 + /(/~-2(l + e)), (5) 
Z 8 = S 2 (l-e) 2 + /(/~-2(l-e)), (6) 

and 

Z 9 = + S 2 Q 2 (1 - e 2 ) 2 + Qp (3 - 1 + e 2 )). (7) 

Intriguingly, the roots of Z-j = correspond to the coordinate singularities associated with 
the event horizon of the KBH, multiplied by 1 + e; for Z 8 = the multiplier is 1 — e. The 
abutment, which lies outside the event horizon, corresponds to a set of orbits for which 
Z 9 = [42J, i.e. 

I 5 + S 2 Q 2 (1 - e 2 ) 2 + Ql 3 (3-I + e 2 ) = 0. (8) 
The solution of equation ([8l is: 



Qx = 2 ^ 2( / 2 _ e2)2 [l (/- e 2 - 3) ± ft (r- e 2 - 3) 2 - 4/(1 - e 2 ) 2 ^ j . (9) 

Where the minus solution must be taken since the plus solution has a singularity at S = 
(unphysical for an SBH) and at e = 1. Further, the value of Q of an evolving orbit cannot 
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exceed Qx] equation (J2J) would yield a complex result for X±. Hence L z and E would possess 
unphysical values. 

By performing an expansion in terms of S 2 one obtains: 

I 2 



Q 



X 



I - e 2 - 3 



Ul-e 2 ) 2 ^ 2 (l-e 2 ) 4 S 4 



I - e 2 - 3 



l-e 2 



5 • • • 



(10) 



Therefore Qx = O \S°) and the j th term of Qx = O ( ) . The expansion of Qx in terms 



of / can be derived from equation (10) once it has been determined to which power of S one 



wishes to work. This result, and its derivatives with respect to I and e, are presented in 



Appendix |A 3| for use in our analysis in section |IV A 



We return to equation (|2]). The terms under the square root can be excluded since Zg = 0. 
Substitution of Qx into the remaining part of the equation yields: 



X% = S 2 



+ 2 



III 2 -4Z-4e 2 + 4 



(11) 



/ - e 2 - 3 , 
2- 10e 2 + e 2 [-3/ + ^ (1 - e 2 ) 2 S 2 



l-e 2 



+ 



6/ 2 + (8e 2 -16)/ + 9-74e 2 + e 4 ) (l-e 2 ) 4 £ 4 



I l-e 2 



From equation ( 11 ) one finds that X± = O I S , I I . Further analysis yields the result: 



X = ±S 



\ 



l(P-4l-Ae 2 + 4\ 



l-e 2 -3 



(12) 



x 1 + 



+ 



l-e 2 ) 2 5 2 



/-e 2 -3) 2 /f[ 2 -4[-4e 2 + 4 



Pi 



'l-e 2 TS 4 



2 ( l-e 2 -3] I 2 [I 2 -Al -Ae 2 + 4 



2^2 
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where 



P 1 = P - (3 - e 2 ) Z + 2 - 10 e 2 , 

P 2 = 5 Z 4 - (34 - 6 e 2 ) Z~ 3 + (84 - 104 e 2 ) Z~ 2 

- (88 - 328 e 2 + 16 e 4 ) Z~+ 32 - 292 e 2 + 200 e 4 - 4 e 6 . 



Orbital energy, E 

The formula for orbital energy, E, for inclined elliptical orbits (see equation (44) in [42]) 
is presented here in a form that more clearly shows that E = O ( S° ) : 



E 



i 



1 - (1 - e 2 



Z~ 3 - Q [I- S 2 ){1- e 2 ) - IX 2 (1 - e 2 



(13) 



further, substitution of Qx ( ) and X\ ( S ) into equation (|13|) yields an expression for E, 



which can be used directly in our analysis, or in the following form: 



S 2 1 - E' 



;i-e 2 ) {t-A)S 2 (1-e 2 ) 2 e 2 l + l-Qe 2 -2)S 4 



l(l-e 2 -3) J2(l_ e 2_ 3 
'1 - e 2 ) 4 (-3 - 30 e 2 + 4 e 2 l + [ 2 + e 4 ") S 6 



P I - e 2 - 3 



(14) 



5. Orbital angle of inclination, i 



The exact formula for i was derived in Boyer-Lindquist coordinates (BL coordinates) and 
found to be: 



Q + L 2 Z + S 2 1 - E 



sin 2 (t) 



Q + L 2 + S 2 (l-E 2 )) -AQS 2 (l-E 2 



2S 2 1 - E 2 



(15) 



which suggests an approximate expansion in the cases of small S or for I — > oo (for which 
E — > 1) |12]. In particular, near-equatorial orbits can also be approximated by such an 
expansion since Q ~ 0. But we are studying near-polar orbits, for which Q > 12; so it is 
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advantageous to exploit the fact that L z = and convert equation (15) to an alternative 
form which can be expanded as a series (not strictly in powers of S 2 ) to obtain: 



cos 2 (<,)=- 



Li ( 1 - E 2 )S 2 



Q + I%-S*(l-&) ( Q + l 2 z -~S 2 (l-E 2 
2Ll (l-E 2 ) 2 > 



+ 



Q + L 2 - S 2 1 - E 2 



(16) 



It is essential to establish the lowest order of S for each term of equation (16); the results 
in equations (10) and (11), and equation (14) to O (S 2 ) will help. 



It was found that X± = — yX± in the vicinity of the abutment (see section 3.5 in 
Komorowski et al. [42j); therefore, 



L K = -JXl + SE. 



(17) 



Each of the expressions in equations (12) and (13), when expanded as a power series in 



I l , will have a leading factor of S and unity, respectively. In evaluating equation (17), the 
leading terms subtract out; therefore, we find that L 2 = O ( S 2 , l~ 2 ) . The inverse dependence 



of L 2 on I is consistent with the physical meaning of L z for orbits on the abutment. Further, 



equation (10) indicates that Qx = O ^S* ,^; therefore, the first term in equation (16) is 
O {s 2 ^ ) an d the second term, O {s^j > with each term containing higher order terms of S 
in increments of 4. 



Taking the square root of both sides of equation (16) yields 



term 



term 



term 



cos (t) = O [S + O [S b )+ O IS 



(18) 



with higher order terms of odd power of S. The second term in equation (16) will contribute 



to equation (18) a factor O \S 5 ); therefore, to derive an analytical formula for (dt/dl 



valid to O (S 3 ) (see equation (1)) it is sufficient to use the first term of equation (16). If 



we choose to work in stronger gravitational fields, for which terms of greater order in S are 



required, then the second and possibly higher order terms in equation (16) would be used 



But we wish to work with terms that contain S and S 3 , to the exclusion of those with S 5 , 



so we shall restrict our analysis to the first term of the series in equation (16); after taking 



the square root, it can be simplified to yield: 



cos (l) 



1^ 
2 Q 



S 2 1 - E> 



Q 



(19) 



Given x = cos (<.), one may calculate i to O ( S 3 ) by using the MacLaurin series for arccos(x) 
to O fx 3 ). 



C. Analytical Formula for i [e,l ) on the Abutment 



We shall now evaluate equation (19) analytically by working with the constituent terms 



as series expansions in 5*, the coefficients of which are expressed in terms of e and /; the 
result to third order in S is our target. An apercu of the method by which the expression 



in equation (19) is treated appears in Appendix B 



1. First- order in S 



To perform our calculation of i to O {^j (see Appendix 



B 1 ) it is sufficient to use: 



(i) _ 7T _ L z 
' = 2 v/Q^ 



where 



(i) 

L z _ _~ s j e2 + 3) (J_ + (1 + e 2 ) + (3 + 2e 2 + e 4 ) (!) -f 5r + 5r'-r'- 



VQx K V/ 3/2 ' l 5/2 I 7 12 l 9/2 

and the number in parenthesis indicates the order in S of the term below it. 



(20) 



(21) 



2. Third- order in S 



Our third-order equations are more complicated. Consider the third-order equation for 



/ : 



(3) - 1 _ _ I 3 

2 6 



(22) 
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where 



in which 



(3) 



x 



V 



X 



L 3 ) \ 



VQ 



X 



J 



(2) \_ 

S 2 ( 1 - E 2 



Q 



x 



(23) 



and 



(3) 



VQ 



X 



(1) 

W 2 + 2 



A 



-5 3 (l-e 2 ) : 



J9/2 



(2) 

/i i 



X 



I 2 I 3 



(see Appendix 



(24) 



(25) 



B2). We evaluate equation (22) to obtain the final result, of [S 3 



(3) 
I 



S 3 (-8-13e 2 -2e 4 + 5/3e 6 ) + (e 2 + 3) (9 + 5 e 2 + 5 e 4 + e 6 ) S 



f-9/2 



+ 



1/2 (1 - e 2 ) (5 - e 2 ) S 3 + (e 2 + 3) (3 + 2e 2 + e 4 ) 5 



f-7/2 



+ S (3 + e 2 ) (1 + e 2 ) T 5 / 2 + 5 (3 + e 2 ) / 



-3/2 _ 



7T 



(26) 



D. Derivatives of Me, Z) on the Abutment 



By taking the partial derivative of i with respect to I (using equation (26)) one obtains: 
dt\ 3 



01 



3 (e 2 + 3) (9 + 5e 2 + 5e 4 + e 6 ) S- (24 + 39 e 2 + 6 e 4 - 5 e 6 ) S 3 



pil/2 



r e 2 + 3 ) ( 3 + 2 e 2 + e 4 ) 5 + \ (l - e 2 ) (5 - e 2 ) S 3 



7 
2 

h - (3 + e 2 ) (1 + e 2 )Sr 7 / 2 - ^ (3 + e 2 )5r 5 / 2 . 



f-9/2 



(27) 



The partial derivative of i with respect to e can also be derived: 

=2e(4(6 + 10e 2 + 6e 4 + e 6 )5- (13 - 5 e 4 + 4 e 2 ) S 3 ) T 9 / 2 
+ 2e ((9 + 10e 2 + 3e 4 ) 5 - (3 - e 2 ) S 3 ) T 7 / 2 

+ 4e (2 + e 2 ) £r 5/2 + 2eSr 3/2 . (28) 
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The formula in equation (27), when evaluated at e = 0, matches the numerical result in 



equation (jlj) for all of the terms with the exception of —122.7 SI 11//2 , which differs slightly 
from the analytical result of —243/2 Sl~ ll l 2 . 



E. Directional Derivatives in the I — e Plane 



Consider the constant of motion, Q, and the corresponding quantity, i, in the I — e plane; 
by using the concept of the directional derivative for two variables, one may represent dQ/dt 
by the equation: 

dQ dQdl dQ de 

~dt ~ ~dldi ~de~dt' ^ ' 

and in a similar manner we may define, 



di dt dl di de 
dt dl dt de dt ' 



(30) 



where the terms dl/dt and de/dt denote the evolution of / and e to arbitrary order. We have 



the benefit of knowing the analytical expressions dQ/dl (see equation (A10)) and dQ/de 



(see equation (All)) at the abutment, which we can derive to the required order. 

A weak- field solution for di/dt, in terms of / and e, has been derived and reported in the 
literature (see equation (15a) in |17|): 



dt mS r ii, 2 ,2 , , 7244 252 
r r^(l - e 2 ) 2 sin (a) ( — + 



19 



e 2 + — e 4 - cos (2^ ) 8e 2 + —e 



26 



(31) 



dt M 2 ' ' ' v / V 15 ' 5 " 2" r», y- ■ 5 

where the term cos(2^ ), in which if) represents the orientation of the elliptical orbit in 
the orbital plane, typically averages to zero with the possible exception where the orbit has 
a large value of e < 1 |17j . More recently, a solution for di/dt to higher order in l" 1 (we 
present the weak-field portion here) was derived by Flanagan and Hinderer [30]: 



di mS r_u 
di = Jl 2 



r^(l-e 2 )^sin(i) 



266 184 
15" " T~ 



151 
~20 



+ — e 2 + ^e 4 + cos (2i) 



22 62 
15 ~ 



39 
20 ' 



(32) 



in which they confirmed a weak-field correspondence to equation (31). In addition to the 



sin (i) contribution found in both equation (31 ) and (32), there is a cos (2t) term in the latter 
expression. 

The trigonometric quantities, sin (a) and cos (a), do not occur in our expressions for i 
and its derivatives at the abutment. But such trigonometric terms are found, usually in 
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a product with S, in the general evolution equations (i.e. dl/dt, de/dt, dQ/dt, and di/dt) 



published in the literature |24"| 1291 l3Tj . One may use equations (20) and (21) to derive 



approximations of sin (i) and cos (t) suitable for working in the leading order of S. Further, 



we may use the approximation of cos (2t) to corroborate the conclusion that equation (32) 



is the same as (31) in the weak-field regime. These trigonometric approximations are only 



valid on the abutment; thus, if it is necessary to perform a differentiation of a trigonometric 



term (as in equation (A14)), then the differentiation must be performed before making the 
approximation. Such limitations notwithstanding, the trigonometric approximations are of 
value to us investigators since they afford us a systematic method for their treatment. 

III. CORRECTION OF di/dl AND di/de FOR SECOND-ORDER EFFECTS 
A. Introduction 

For circular orbits, Komorowski et al. found that the numerical estimate of 

(di/dl) in the weak-field regime deviates from the di/dl results reported in the liter- 

V / min 

ature (see Flanagan and Hinderer [30], and Hughes |21j). Consider the quotient of the 
formulae presented in equation (3.9) of Hughes (21] where i = tt/2: 

W_ = ^ = _61^_|_ (33) 

Rweak 01 48 

Because —61/48 > —4.5 in the weak-field regime, X\ =>• X 2 is the pertinent mode; and 
the best information one can obtain from idi/dl ) is the specification of the lower limit 

V / min 

of di/dl for all I > Ilso, abutment ( see |12] for more details about the abutment and its 
relationship to the last stable orbit (LSO)). Therefore the second-order (i.e. d 2 Q pat h/dl 2 ) 
behaviour at the point of tangential intersection of Qx and Q pa th must be considered. In 
section [IT] the numerical results have been verified by analytical derivation of the formula 
for (di/dl) to O ( S 3 ) . It remains for us to extend this analysis to include second-order 

V / min V / 

effects on elliptical orbits; to this end, we shall discuss how to incorporate second-order effects 



into Qx, and the resultant change to the formula for X± (see equation (j2J)). Equation (20) 
is sufficient in treating Xl, and then ultimately i, to the leading order in S. 
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B. Second-order Effects in Q pa th 



1. Circular orbits 



Let us begin our treatment in the Q — I plane with the value of e held constant at zero. 



The form of Qx is represented by the series in equation (10); and because Q pa th intersects 
Qx tangentially at a single point (7 ) (contact of the first order), we surmise: 



Q 



path 



Q 



X 



l=lo 



(34) 



and 



dQpath 



dl 



dQ: 



l=lo 



dl 



(35) 



l=lo 



But the abutment can only offer an upper bound on the second derivative of Q pa th, i-e. 

d 2 Qpath 



dl 2 



< 



d 2 Q x 



l = lo 



dp 



(36) 



l = lo 



To perform an analytical treatment of the second derivative of Q pat h, we define an ansatz 



Q 



path 



Q 



A 5 



X 



I (I) 



where 



f(lo 




(37) 



(38) 



and 



a > 0, 



X — l-L 



(39) 



where p and q shall be determined by requiring that the weak-field solution be of the form, 
SI" 2 (see equations (27) and (33)). The adjustment represented by equation (37) is based 



on the Taylor expansion of a function; the function / yl j represents a second derivative of 
a primitive, p(e, I), with respect to /, which is evaluated at l . One must not confuse the 



concepts of the abutment and Taylor series; equation (37) is not intended to be a Taylor 



series representation of Q pa th- We have taken the analytical formula for Qx and incorporated 
a term, which is designed to adjust the second derivative of Q pa th so that it makes contact 
with Q x tangentially at a prescribed point, l Q . If I = l Q , then the adjustment to Q pa th an d 
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dQpath/dl is zero; and the value of d 2 Q pa th/dl 2 is reduced by / u 1. Equation (37) can be 



applied to the analytical development of di/dl. We shall call this reduction of the second 
derivative the reductive ansatz circular. 

Let us consider how the reductive ansatz circular affects equation ([2]), with attention 
given to equations ([7]) and ([8]); although we begin with an analysis of circular orbits, e has 
been retained in these equations for later use. Evaluate 



Z 9 (Qpath) = Z 9 (Q X - $) 

= f + S 2 (l-e 2 ) 2 Q 2 x - (1-3 " 2 ^ 3 
+ $ (2S 2 (1 - e 2 ) 2 Q x + $S 2 (1 - e 2 ) z + ll- 3 - e 2 ) I 3 



<-\TQx 
, 2 



(40) 



where 



for which the quantity in square brackets in equation (40) is equal to zero (viz. equation 



abutment 



; therefore, the use of this reductive ansatz has assured 



8J)) for all values of I > Ilso, 
us of an effective means to simplify the expressions. The terms that remain share a common 
factor, A 2 , which will appear as ±A when taken outside of the square root in equation 
We shall limit our analysis to O (s^j (the S 2 terms will affect terms of higher order in I" 1 
in the series in equations (27) and (28)); therefore, the product of Z 7 (equation (|5|), Z 8 



(equation Q), and Zg (equation (40)) simplifies to: 



Zy/i — ZiZ^Zt 



'5 



$f [I - 2 (1 + e) )[l - 2 (1 - e) ) [I - 3 - e 



(41) 



where we evaluate Z^ at the point of intersection on the abutment by setting l Q = I (i.e. 
A = 0). We take the square root of Z/^, and a term, ±Aa/2/2, emerges. The choice of sign 
is determined by the mode at the abutment. 

Mathematically speaking there are two modes at the abutment: the fast mode 



X 2 => xl, 



(42) 



and the slow mode 



Xl 



X 2 



(43) 
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In section 5 of [42] it was established that orbits that evolve on a path towards the abutment 
(during which I > l and A > 0) are governed by Xi (see equation (2)) and after making 
contact with the abutment at I = l the orbits are then governed by X 2 _ (for which I < l 
and A < 0) (see figure [I]). Thus by choosing the positive sign for ±A the equation remains 
consistent with the dominance of the slow mode. If one were to perform an analysis for the 
fast mode then —A would be used instead. 

An examination of equations (|3]), (HJ), and (10) reveals that 

Z 5 + Z 6 Q X = 0(S 2 J 4 ), (44) 



from which one may infer 



2SJZ^ = 0[S 



V2XS x S q/2 = O (S 2 ) ; 



(45) 



(46) 



therefore, q = 2 in the reductive ansatz (see equation (38)). The value of p can be derived 
by considering the order of L z in I. We find (viz. equation (17)) that L z = O ), which 



must not be changed by the reductive ansatze. And the leading term, S, in the expression 



for X (see equation (12)) must remain. Given the order of I in equation (44), one must work 

(47) 



with the next lower order, i.e., 



2SJZ 



Oil 



± ^?A/ p / 2 / 4 
2 



0(1 



(48) 



Given A = O ylj , we conclude that p = —4. In our reductive ansatze, we have found the 
values of p and q that ensure the second-order effect does not change the form of dt/dl in 
the weak-field regime. 



2. Elliptical orbits 



The general formulation of the reductive ansatz elliptical can be derived by starting with 



a Taylor series for two variables (see Appendix Al). Because we concern ourselves with 
second-order effects, we shall use the following operator: 



1 

2! 



x d d 
A— ~ + e^ 
dl de 



(49) 
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where A = (/ — l a ) and e = (e — e D ), and where the ordered pair ( e , l a J specifies the location 
of the contact of first order between Q pa th and Qx (see figure [TJ . 
One may define the reductive ansatz elliptical, i.e. 

2 



Qpath — Qx 

= Qx 



l 

2 
A 2 



A^ + e — 
9/ <9e 



<9 2 

— p(e, /) 
dP 



p(e,l) 
- Ae 



J 6 — 605/ — /o 

<9 2 



6 — 6o)^ — 



a_ 2 

<9e 



(50) 



where we conjecture the existence of a primitive function, p(e,l). Equation (36) will also be 
applied to the case of elliptical orbits where e is close to zero. 

The expression for Q pa th is best regarded as a parameterized curve, and to make such a 



treatment in equation (|50|), one may factor out the A, to obtain 
A 2 



Qpath — Qx 

= Qx 



2 

A 2 



d 2 

-^p(e, I) 
dl 2 



+ 2 



d 2 



dlde 



p(e,l) 



+ 



de 2 



p(e,l) 



9 e , l 



(51) 



for which we have the benefit of knowing the limiting form of e/A (= de/dl) to arbitrary order 
in l~ x . Thus it is possible to parameterize Q pa th in terms of A. We can use the expression, 

^a t (e ) (52) 



g (e , / ) =S 2 J2y t 



;i+4 
i=0 l o 



and this will form the basis of the reductive ansatz elliptical. 



C. Application of the Reductive Ansatze to the Analytical Derivation of dijdl and 



di/de 



The reductive ansatze (equations (37 39) and (equation (|5Tj)) constitute a reduction of the 
second derivative of Qx to more realistically model the behaviour of Q pa th at the abutment 
and perform a methodical treatment of this reduction in the analytical calculation of dijdl 
and di/de. 



The procedure outlined in Appendix A 2 yields the following formula for dijdl: 

d I {r.l.X.S) d 



dl 



A=0 



dl 



i \ e,l, A, S 
~d\ 



dX 
dl 



(53) 



A=0 
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where 



ax 



1 and A 



l=lo 



but the result for dt/de is simpler, 



dt 
de 







i[e,l,\,S 



A=0 



de 



de 



(54) 



Equation (54) and the first term in equation (53) yield the formulae that describe the 
evolution of i for a Qpath along the abutment (i.e. (<9t/<9e) min and (di/dl) min ). The second 



term of equation (53) incorporates second-order effects, and thus describes the physically 
more realistic situation in which Q pa th intersects the abutment tangentially at a single point. 
Because one takes the first derivative with respect to A, the second and higher powers of A 
will vanish when setting A = 0. But as we shall presently see, the second-order effects of the 
reductive ansatze remain. 

We choose to work with the symbols e and I rather than e a and l a , given that e and 



I can be used to represent an arbitrary point on the abutment (see Appendix A 2). The 
reductive ansatz circular (see equations (37 39)) is applied at the abutment with p = —4, 
q = 2, and n = 2 (while retaining the two terms of leading order in S at the conclusion of 
the calculation) with 



A)(o) r 5/2 



(55) 



where 



a i \ / — -— , . ai(e) 

^o(e) = + — V«o(e), and Ay (e) = — — - 



To apply this method to elliptical orbits, we will be required to use g ( e, / ) . To calculate 
that function, the common primitive p(e, I) is needed. 



D. Analytical Derivation of the Common Primitive, p(e, l) 

Now that the values of the parameters, p = —4 and q = 2, have been found, it is possible 
to derive the formula for p(e,l). Consider the reductive ansatz circular: 

n 

/(0=« 2 E^- o*> 

i=0 1 
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We conjecture a more general form of / (t\ that includes e: 



n 

/(«.9 = s, £ 

i=0 



di (1 + hie 2 



li+4 



(57) 



Performing the first integration over I yields: 



F ej 



f[e,l)dl 



-S 2 



E 

i=0 



1 a,; (1 + hie 1 
i + 3 fi+3 



x(e) 



(58) 



The second integration over / yields an expression for the common primitive: 



p{e,l) = I F(e,l)dl 



S 2 



n 

E 

,i=0 



di (1 + 6je 2 



+ 2) (z + 3) £+2 



+ x(e)[+C(e) 



(59) 



The constants of integration, x(e) I and ( (e), can each be set to zero since we require 
lim p(e, I) = 0. Now that the formula for p(e, I) is known, it is possible to obtain g(e,l), 
which is required by the reductive ansatz elliptical. 



IV. THE TREATMENT OF dQ/dt AND dijdt ON THE ABUTMENT 



A. The dQ/dt Evolution Equations 



Komorowski et al. [12] investigated the consistency of dQ/dt with the evolution equation 
dl/dt, for circular orbits at the abutment (l g2 7r/2) by performing a preliminary numerical 
analysis for values of / = {7.0, 100.0} and KBH spin S = {0.05,0.95} (see section 5.2.1 of 
|42j). The published values of dl/dt |21j . which we used in our investigation, were calculated 
for i ~ 7r/3, and the difference of this value of i from that at the abutment contributed to 
some inaccuracy in the analysis [12]. In this work, the derivation of analytical formulae for i 
and its derivatives, as well as the use of the directional derivative to determine dQ/dt, now 
allow one to perform a more complete treatment for elliptical orbits. 



Let us consider the directional derivative in equation (29) as a means of deriving dQ/dt 
at the abutment. We have demonstrated that the second-order effects are not seen when 



calculating the first derivatives of Q pat h (i-e. dQ pat h/dl and dQ pat h/de, see equation (35) or 
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(51)); therefore, we may use dQx/dl (equation (A10)) and dQx/de (equation (All)) when 



working with equation (29). 



The form of dQ/dt (equation (A. 3) in [28j (after equation (56) in [27])), which was used 
in |42j to test dQ/dt will be revisited in this work: 

_ / _ lS 2 (3 + e 2 f\ 64 m 2 

2PN \ ^ 



dQ 
dt 



(i C 2\V^VQ 

5 M v ; [7/2 



011 , 012 

Oq = h 7T~ 

11 £3/2 



, §2 9w (3 + e 2 ) | 45g 4 (3 + e 



I 2 

2^2 



(60) 

P 8 Z 5 J 

But it is preferable that the formula for dQ/dt (and for di/dt) that we test be accompanied, 
in the same work, by their associated expressions for dl/dt and de/dt; and fortunately a 
paper by Ganz et al. |29| provides such information, which we shall use in our analysis. In 
particular, we will use equation (4.3) in [29]: the evolution equation for I, 



dt 



X 



fl . #12 . f% — fiS 2 f 7 
Oq + 71- 1 = 71- 

I Pi 2 I 2 I 5 / 2 



^ f 2 (3 + e 2 )S 2 f 6 (3 + e 2 )S 2 | / 5 (3 + e 2 ) 2 S 4 



P 



P 



P 



(61) 



which is (excluding the common factor, / 3 ) to O (^l 5 / 2 j in [29] : and the evolution equation 
for e, 



de _ _304 / _mx ~ 4 , 
dt ~ 15 \MV [ > 



h 2 h± h 5 + h 6 S 2 hq 



h\ —-\-7i- 

l P/2 P 



— 71 — 

15/2 



^ h 3 (3 + e 2 )S 2 h 8 (3 + e 2 ) S 2 | h 7 (3 + e 2 ) 2 S 4 



P 



P 



P 



(62) 



which is (excluding the common factor, t 4 ) also to O (r 5/2 ) in [29 



20 



The evolution equation of Q 

64 / m 



dQ 
dt 



2.5PN 



04 / m \ ~ , , o\ ; , 



S 2 (3 + e 



X [fl-9 



rfl , #12 4 - ^5' 2 ^7 

h 7T~ ~ 7T~ 

I 1 /3/2 /2 J5/2 



S 2 (3 + e 2 ) rf 2 S 2 (3 + e 2 ) 4 S 4 (3 + e 2 ) 2 d, 



/ 3 



+ 



(63) 



which corresponds to equation (4.1) in |29j . was of O (^l 5 / 2 j (excluding the common fac- 
tor, l~ 3 ). These formulae (equations (61), (62), and (63)) have been converted from the 



variables used in |29j to our variables (see Appendices A 4 and A 5). Because some of the 



original coefficients contained cos (l), which we have replaced with the approximation based 
on equation (21), there are new terms, which correspond to Z -3 , Z~ 4 , and l~ 5 , in each of 



equations (61), (62), and (63). The original expressions did not include terms with these 



powers of I, so we cannot use the new terms to extend the accuracy of our analysis beyond 
that of the original expressions in Ganz et al. [21] • Further, these evolution equations are 
O (e 2 ), hence the final results must also be used up to the second power of e. 



We assume that the evolution of the orbit, dl/dt and de/dt, is described by equations (61 ) 



and (62) (from Ganz et al. [29J); the result of evaluating equation (29) is compiled in table 
[T] (first column). The second column contains the result (dQ/dt) 25PN , derived by Ganz et 



al. (see equation (63)), evaluated on the abutment. Similarly, the third column contains the 

also evaluated on the abutment. 



formula for (dQ/dt) 2PN based on equation (A. 3) in 

Although the terms in table [i] are O (e 4 ) (N.B. Those given to O (e 2 ) were found to have 
no terms of higher order in e 2 ), any final result one might calculate must be of O (e 2 ). To 
that order, our calculated results at the abutment agree with those of Ganz et al. [29] up 
to /~ 5//2 . There is also agreement with Barausse, Hughes, and Rezzolla [28J up to l~ 2 with 
the exception of the coefficient for the 5* 2 term in the third column (marked with §), which 
differs from the other two results; the expanded equation in Barausse, Hughes, and Rezzolla 
differed from that of Ganz for that order of I. 

There are two reasons for reporting these results to O (e 4 ): first, we wish to demonstrate 
that the confirmation of Ganz's calculations is not to be dismissed as a fortuitous triviality. 
The method of calculation of dQ/dt at the abutment differs fundamentally from that used 
by Ganz et al. to derive their results, as would be required of a good consistency condition; 
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Table I: The coefficients of —5/64 M 2 /m — e 2 )~ 3 / 2 dQ/dt. The first column contains our 
calculated results at the abutment. The second column contains the results of Ganz et al. [29J, and 
the third column contains the results of Barausse, Hughes, and Rezzolla [28J. The trigonometric 
functions in both of these sets were evaluated on the abutment. Note: although the terms are 
reported to 0(e 4 ) they are only accurate to 0(e 2 ). 



Results at the Abutment 
(see equation (29)) 



(equation 4.1 in |29| 
J) 



(see equation 



68 



(equation (A3) in |25] ) 



(see equation (60)) 
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second, the differing coefficient values for the e 4 terms demonstrate that one must not per- 
form calculations on the abutment for highly eccentric orbits. While the abutment equations 
{dQx/dl, dQx/de, and t (l, e, Sj) are exact in terms of e 2 , one remains limited by the order 
of e used in the radiation back-reaction model being tested. 

Because the expressions for dQx/dl and dQx/de can be derived to arbitrary order in Z -1 , 
and the coefficients for each power are exact finite series in e 2 , it is worthwhile to consider 
using the abutment to improve the order of e 2 of the evolution equations in the weaker 
field regime by allowing other theoreticians to perform a test of their own, improved back- 
reaction models. Since the abutment extends down to the LSO, one might also explore the 
development and testing of evolution equations in the strong-field regime, given that on the 
abutment the trigonometric contributions of sin (t) and cos (l) can be expressed as functions 
of e, I, and S. But one must also be mindful of the assumptions made at the outset of this 
exercise, in particular, the assumption that the secondary object can be approximated as a 
test-particle of infinitesimal mass, and the use of adiabatically evolving orbits. 
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B. The Second-order Calculation of di/dt for the Leading Order of S (weak-field 
regime) 

Now that p(e,l) is known we can calculate g (e,t\; but let us first derive dt/dl using 



equations (61) and (62). We find the following: 



de 19 n-me 2 3215 33373 - 2 



dl 12 



304 _|_ 3192 102144 



I 



/ 2 



expressed to O [I 2 ). From equation (51), one derives: 



glej) = S 



a [e) a x (ej 



where 



and 



I 4 I 5 
119 

a (e) = oo (0) ( 1 - ^b e 2 



, J ( 3215 , 143 \ 2 \ 

ai ( e ) = ai (0) 1 + b &i e 2 , 

v ; \ V72576 864 J J ' 

which can be used to calculate dt/dl under the reductive ansatz elliptical, 

dL -S { I (3 + e 2 ) (1 + e 2 ) + \ (3 + e 2 ) A (e) - i ^ (e) } T 7 / 2 
-^{^(3 + e 2 )-A ( e )}r 5 / 2 , 



9/ 



where 



and 



A (e) = A (0) ( 1 - l ~b Q e 2 



A 1 (e) 



ai (e) 
A (el 



(64) 



(65) 



(66) 



Now that we have developed a formula for dt/dl (equation (66)) that incorporates the 



reductive ansatz elliptical, and we have found that dt/de (equation(28)) is unaffected by the 



reductive ansatz elliptical, the expression for di/dt can be obtained from equation (30) to 
the leading order in S with coefficients of O (e 2 ): 



£™r 4 (l-e 2 )^ u * 
M 2 v ; V/ 3 



di ~ m r 
dt 



[3/2 [5/2/' 



(67) 
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where 



32 4 
U x = — (9 - 2A (e)) + — (-109 + 42A (e)) e 2 
15 15 

^3 = 77^ ( 1647 - 2494A o ( e ) + 168A i ( e )) 
105 

+ — (682 - 2978A (e) + U7A X (e)) e 2 . 
105 



Equation (A15) can be expanded and expressed to leading order in S to yield: 



di 244 ~ m r 



dt 15 M 2 1 ' I [3/2 [5/2 r 



(68) 



By equating the terms in equations Ip7\) and (68) (i.e. XJ\ = 244^/15 and U3 = 244m 3 /15) 



one first solves for A (0) and A\ (0) for a circular orbit by setting e = 0: 

155 



A) (0) 



18 



and 



A (0) 



279289 
4032 



By substituting these values into equation (67), we obtain 



di_ 244 5_m_r_4 

dt~ 15 M 2 [L 6) 



(l + ( 



18115/) 

52704 488 



213) e 2) 



(10461 1 / 
V 1708 ^ V 



/3/2 

79869 _ 39938327 
54656 17708544 



6l + 



5621763839 
1487517696 



bo) e 2 ) 



15/2 



(69) 



When evaluated at e = 0, the expression in equation (69) matches the results reported in 



the literature (equations (31), (32), and (68)) for near-polar orbits. For near-circular orbits, 



values of 6q an d b\ can be found for which the coefficients of the I 3//2 and I 5 / 2 terms in 



equation (69) match their theoretical counterparts in equation (68). 



C. The Independence of the Abutment of Radiation Back-reaction Models 

Let us clarify the meaning of our statement that the abutment model is independent of 
any specific radiation back-reaction model. The expression for the abutment, Qx (equation 
rt9J)), is determined by the characteristics of the Kerr spacetime of the primary object in 
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which the secondary object (i.e. test-particle) orbits. The analytical expressions for dl/dt 
and de/dt describe the effects of radiation back- reaction on the values of I and e of the orbit, 
and they serve as inputs to our abutment model in two ways: first, through the quotient 



e/A = de/dl (equation (51)); and second, through the directional derivatives in equations 



(|29]) and (go]). 

The mechanics of the abutment remain consistent, the details of the radiation back- 
reaction model notwithstanding. The results of either directional derivative are outputs of 
the abutment model that describe the effect of the radiation back-reaction on the listing of 
the test-particle orbit. 



V. CONCLUSIONS 

For inclined test-particle orbits around a black hole, two solutions for X 2 (where X = 
L Z — SE) can be derived: X 2 _ and X?. Given a Schwarzschild black hole (SBH), X^_ = X 2 ^ on 
any polar orbit, where X 2 corresponds to prograde orbits and X\ corresponds to retrograde 
orbits. For a Kerr black hole (KBH) the orbits on which X^_ = X 2 _ are not polar, but 
near-polar and retrograde. Such orbits comprise the abutment at which the value of the 
Carter constant (Q) is a maximum for given values of latus rectum (I) and eccentricity (e). 

In this work we derived an analytical formula for the value of orbital inclination, l, of an 
elliptical orbit on the abutment. By performing the partial differentiation of i with respect 
to I, we were able to confirm the numerical result for dt/dl reported in Komorowski et al. 
|42| for circular orbits, and we were able to extend the formula to include dt/dl for elliptical 
orbits. A result for dijde was also obtained for elliptical orbits. Further, it allowed one 
to redefine, in terms of e, /, and S, any trigonometric function that might be found in an 
evolution equation to be tested at the abutment. 

Evolving orbits in Kerr spacetime are not constrained to follow the abutment. Instead, 
the value of Q will follow Q pa th, which intersects the abutment tangentially at an arbitrary 
point of contact of the first order. This behaviour is assured because the value of Q pa th cannot 
exceed that of Q on the abutment; to do so would make X± complex and thus unphysical. 
For circular orbits, we modelled the second-order behaviour reported in [42J by introducing 
a bounded function / ^e, (also in terms of S) to reduce the value of d 2 Q pat h/dl 2 while 
leaving Q pa th and dQ pat h/dl equal to their corresponding values (Qx and dQx/dl) on the 
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abutment. This approach was then applied to elliptical orbits, and a new bounded function 
g (e, l\ , which depends upon de / dl, was used to reduce the value of d 2 Qx / dl 2 to d 2 Q pat h/ dl 2 . 
It was discovered that the value of d 2 Q pat h/de 2 remained unchanged by the reductive ansatz 
elliptical. 

The consistency of published evolution equations, dQ/dt, dl/dt, and de/dt, was tested by 
using dl/dt and de/dt to generate an expression for dQ/dt at the abutment. In general, the 
calculation of dQ/dt is more difficult to perform than that of dl/dt and de/dt [19|; hence, 
the abutment provides a useful mechanism for testing the validity of radiation back-reaction 
models. Indeed, the evolution equations reported by Ganz et al. [29| were confirmed to 
their 2.5PN order. The abutment provides a consistency condition that is limited to near- 
polar retrograde orbits; and yet, some back-reaction models have been found to exhibit 
pathological behaviour for polar orbits (Gair and Glampedakis [27J). Another consistency 
condition is already known for orbits of any i (see equations (13) and (14) in [27]), but 
it applies to circular orbits (i.e. in the limit e — > 0); the abutment is valid for orbits of 
arbitrary eccentricity, depending on the accuracy of the back-reaction model used. 

This method promises to be a useful tool for confirming the accuracy of evolution equa- 
tions to greater order in e and l~ x . Further work might entail the development of a more 
precise mathematical treatment of the ansatze in relation to the underlying physical concepts 
of the radiation back-reaction process and its effect on the listing behaviour of orbits near 
the abutment. It would also be intriguing to investigate the consistency condition reported 
by Gair and Glampedakis |27j, on the abutment. 
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Appendix A: Ancillary Equations 



1. Taylor Series for two Variables 

Refer to Chapter 6 in |44j for a more detailed treatment. Let us consider a locally 
continuous function with two independent variables, f(x,y). We may use an operator 



h— + k— 
dx dy 



(Al) 



to construct a Taylor series of n terms 



/ (x + h,y Q + k) = f (x , y Q ) 
+ 



h Tx + k lfy) nX ' y \ 



+ 



x=x ,y=y 



2! 



h iL + k iry ] f{x ' y) 



... + 



TV. 



r\ rj \ n 



(A2) 



x=x ,y=y 



if the (n + l) th partial derivatives are continuous. In this paper, we are concerned only with 
the second derivative. 

2. Treatment of the Taylor Series Under Partial Differentiation 



Given the term: 

A = hf (x )g(x), 

where h = (x — x Q ). We can calculate the partial derivative of A with respect to x, 

dA d 



(A3) 



dx dx 



(hf (x )g(x)) 
d 

f(xo) (hg 0)) 

f d d 
f {xo) ( h—g (x) + g (x) —h ) , 



(A4) 



and thus demonstrate 



dA 
dx 



f (x ) g (x ) 



(A5) 



27 



Consider a more complicated case where we have a function F(x, h,y, k), where h 
(x — x ) and k = (y — y Q ). Calculate 

dF 



dF 

dx 



If we hold y constant and set k — 0, then 



and 



dy 



y=yo 



dF _ dF dF dh 

dx dx dh dx 



If we hold x constant and set h — 0, then 



(A6) 



(AT) 



dF _ dF dF dk 
dy dy dk dy 



(A8) 



These results will be of use in applying the second-order effects to Q pa th as it makes contact 
with the abutment, Qx- 

3. Treatment of Qx as a Series in I 



The expansion of Qx in terms of S (equation (10)) can be expressed as a series in I: 



Qx = UJ2 



(3 + e 2 ) m , i (i - 1) (3 + e 2 )^ (1 - e 2 )" ~ 



+ 



s 1 



(A9) 



From equation (A9) one can obtain dQx/dl and dQx/de directly: 
dQ x 



and 



di 



dQ 



i=0 



i (3 + e 2 ) l+1 i 2 (i - 1) (3 + e 2 r 2 (1 - e 2 f ~ 



Z* +1 



+ 



x 



de 



2e E 



(z + l)(3 + e 2 y 



z(z-l) 



li+i 



(3 + e 2 r 2 (1 - e 2 



5 2 



(A10) 



-S 2 



»(»-!)(» -2) (3 + e 2 r 3 (l-e 2 ) 2 ^ 

2 r* 



(All) 



28 



The y/Qx will also be required for the treatment of dQ/dt, (see Appendix A4) 

*2 ^9 , K /r. , 2\3 



/- - l(3 + e 2 ) , 3(3 + e ») a , l ((l-e 2 )^ 2 + |(3 + e 
5 (3 + e 2 ) (32 (1 -e 2 fS 2 + 7 (3 + e 



J5/2 



+ 



+ 



128 F/2 
7 (3 + e 2 ) 2 (80 (1 -e 2 ) 2 S 2 + 9 (3 + e 



2^3 



256 



p/2 



+ ... 



(A12) 



4. The 2PN Flux for Q 



Equation (60), an expression for dQ/dt, was derived from equation (A. 3) in [28J (af- 
ter equation (56) in [27]) by substituting the approximations of sin (a) and cos (l) on the 
abutment, where 



99 = 1 + 8 6 ' 9w 

a , 97 - 
#12 = 4 + — e 



7 2 -6 _ 61 91- 461 4 
+ 4 + 64 ' 

44711 302893 
+ 



1247 425 



9u 



913 



9u 



336 + 336 6 ' 
33 95 2 



9072 ' 6048 ' 16 ' 16 

An alternative expression for dQ/dt to 2.5PN order was presented by Ganz et al. (equa- 
tion (4.1) in [29]) in which we have converted their variable, Y ~ cos (t), to — S (3 + e 2 ) Z -3 / 2 



to yield equation (63), where 



743 23 2 , 85 211 „ , 

di = e , a 2 — 1 e "a 

336 42 ' 



129193 84035 
+ 



18144 172 



329 929 

«4 — ~7Z7T + — e 



~ 2 , G?^ 

96 96 ' 5 



53 163 2 , 2553 553 2 , 4159 21229 2 
— H — — e , 06 = -— — T7^ e > "7 ~ 



224 192 



+ 



672 1344 



e . 



5. Evolution Equations for I, e, and t 

The evolution equations for I, e , and 1 are reported by Ganz et al. (see equation (4.3) 



in [29]) to O 1 1 5 / 2 ). We reproduce them in equations (61), (62), and (A15) after having 



converted their original variable, 


v = 


yff, to t 1 / 2 , 


743 55 2 

A = 1 e 

J 336 21 ' 


h = 


133 379 2 
l2" + ^4" e ' 


329 929 2 

U = 1 e , 

J 96 96 


h 


815 477 

= 1 e 

96 32 



34103 526955 



18144 12096 ' 
1451 1043 o 



56 



+ 



96 



29 



4159 48809 2 

It = 1 e 

J7 672 1344 ' 



and where 



121 , , 6849 4509 , , 879 515 2 

hi — 1 H e , /io = 1 e , % = 1 e , 

304 ' 2128 2128 ' 76 76 ' 

985 5969 2 , 286397 2064415 2 , 3179 8925 2 

/u = 1 e , n.5 = 1 e , ha = 1 e , 

152 608 38304 51072 ' 608 1216 

5869 10747 2 , 1903 22373 2 , 87947 4072433 2 

Ai? = 1 e , /is = e , /ig = 1 e . 

608 1216 304 8512 ' 4256 68096 

We use the following equations to convert the form of the equation for list rate in |29j: 

1 - Y 2 ~ sin 2 (i) ; (A13) 



We obtain: 



c?F dcos(i) di 

dt di dt 

, s dt 



— sin U x 



(A14) 



di 244 / m 



dt 15 \MV 



1 



r 4 (i-e 2 ) ! (i-^ 2 (3 + e 2 ) 2 r 3 



Ul S _u 3 S u 2 (3 + e 2 )S 3 
244 / m 



15 (S)^ 1 -^ 

ui5 u 3 S u 2 {3 + e 2 )5 3 



/3/2 J5/2 



/7/2 



1 m (3 + e 2 ) 2 g 3 l ^ 3 (3 + e 2 ) 2 g 3 1 n 2 (3 + e 2 ) 3 ^ 5 ' 

2 [9/2 + 2 /H/2 2 /13/2 



(A15) 



where 



Mi = 1 + 



189 



13 277 



61 " ' U2 244 ' 244" ' 



«3 



10461 83723 
+ 



1708 3416 



Appendix B: Series Expansions of Critical Values in Terms of S and I 



Our conversion of the quantities X (equation (12)), E (equation (13)), and Qx (equation 



(10)) to expansion series in S helped to simplify our analysis by avoiding the use of the 
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much more complicated series expansions in terms of I. Equation (19) can be converted to 
a series: 

(2n+l) 

f - 

J> m ^ +1 . (Bl) 



By choosing the order of S (the value of 2n + 1) in which to work, it becomes easier to derive 
suitable series approximations of these quantities, and their mathematical combinations, in 
terms of I. Since the equations derived during the full analytical treatment are Brobding- 
nagian, and thus preclude detailed presentation in this paper, we shall offer the essential 
highlights of our analysis. 



1. First-order Calculations 



We require the series expansion of the quotient, which appears in equation (20) 



(i) 
I, 



(B2) 



x 



to be expressed in terms of /. To obtain this result we perform a careful manipulation of 



L z (in terms of X and E, viz. equation (17)) and Qx (as a series expansion in S) using 



MacLaurin series. The coefficient of S 1 (i.e. c\) is converted to an expansion in /; we find 
Ci to be: 



1p! - 1p 2 



l-e 2 -3 



where 



01 



i 



l{ l 2 - 4 [ I + e 2 - 1 



Z-e 2 



^2 



\ 



I 2 - 4 / + e 



I Z-e 2 



From equation (|B3|), one obtains the result 
n = - ' + 1 — + — r— + 



(e 2 + 3 ) (—+ (1 + g2) ■ (3 + 2e 2 + e 4 ) , (9 + 5 e 2 + 5 e 4 + e 6 ^ 



V/3/2 



/5/2 



/7/2 



p/2 



(B3) 



(B4) 



which appears in equation (21). 
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2. Third-order Calculations 



The third-order calculations require two additional factors: 



(3) 

I, 



VQ 



and 



(2) 

s 2 f i eA 



X 



Q 



X 



(B5) 



which are used to evaluate x = cos(t) using equation (23). The first factor can be derived 
by converting the coefficient of S 3 , 



:i- e 



2N2 



c 3 = — 



7 2 



l 2 -2l + 2le -16e^ Vi+ ~ 4 + 1 ) \ l ~ 2 ~ 2 e ]i>: 



l-e 2 -3 



l(l-e 2 -3) / 2 -4/ + 4e 2 -4 



(B6) 



in equation (Bl), to a series expansion in I (see equation (24)) and adding the result to the 
first order term: 

. o / 1 1 1 1 4- F> \ 

(B7) 



c 3 



V ^ ^7/2 2 /9/2 y 



see equation (24)). The second factor in equation (B5) is also obtained by working in 



expansions of S, which proceeds by a simpler derivation (see equation (25)). The orbital 



inclination, l, is then obtained by using equation (22) 
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e=e r 



r 



Qpath = Qx _ 



Qpath < Qx 



X 



Qpath < Qx 



l=l Q 



1 



Figure 1: A generic representation of Q pa th (short-dashed line) is depicted in the I — e plane as 
it makes contact of the first order at a point (Z ,e ) on the Qx surface. The direction in which 
the orbit evolves is shown by the arrow. We offer a generic representation to emphasise that any 
Qpath, which is predicted by a radiation back-reaction model, may be tested in this manner. 
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